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Abstract — Optical  supermode  dynamics  are  nowadays  receiving 
an  ever  increasing  attention  in  mode-division  multiplexing  appli¬ 
cations.  In  this  paper,  we  systematically  analyze  and  characterize 
these  modes  in  various  structures  based  on  coupled-mode  theory. 
In  particular,  we  investigate  the  structure  of  supermodes  in  hon¬ 
eycomb  multi-ring  arrangements  of  any  size.  In  addition,  we  show 
that  higher  order  supermodes  in  coupled  few-mode  multi-core 
waveguide  arrays  can  be  strongly  affected  by  angle-dependent  cou¬ 
plings,  leading  to  different  modal  field  profiles.  Analytical  solutions 
are  provided  for  linear,  rectangular,  and  ring  arrays.  Higher  order 
supermodes  are  observed  for  the  first  time  in  a  coupled  few-mode 
three-core  fiber  system  using  the  S2  imaging  method. 

Index  Terms — Microstructure,  optical  waveguides,  optical 
fibers,  optical  fiber  communication,  space  division  multiplexing. 


I.  Introduction 

SUPERMODES  are  eigenmodes  of  composite  structures  in¬ 
volving  coupled  constituent  elements,  each  of  which  also 
supporting  guided  modes  in  isolation  [1],  [2].  As  indicated  in 
a  number  of  studies,  such  multi-core  systems  can  be  effec¬ 
tively  analyzed  using  coupled-mode  theory  (CMT) — which  is 
particularly  effective  when  the  coupling  between  neighboring 
elements  is  relatively  weak.  Over  the  years,  the  properties  of 
supermodes  in  either  linear  [3]  or  ring  arrays  of  coupled  waveg¬ 
uides  have  been  analyzed  using  CMT  methods  [4]— [6] .  Even 
in  the  simplest  possible  configuration  of  two  coupled  channels, 
supermodes  play  an  important  role  given  that  their  interference 
is  the  one  responsible  for  the  energy  exchange  behavior  in  a 
directional  coupler. 

Quite  recently,  supermodes  have  received  renewed  interest 
within  the  context  of  mode-division  multiplexing  (MDM),  a  new 
transmission  method  aimed  at  overcoming  the  capacity  limit  of 
single-mode  fiber  communication  systems.  Since  modes  tend 
to  couple  during  long-distance  fiber  transmission,  unravelling 
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mode  crosstalk  using  multiple-input-multiple-output  (MIMO) 
digital  signal  processing  (DSP)  is  necessary  for  demultiplexing 
MDM  channels  [7] .  In  the  presence  of  modal  group  dispersion, 
the  computational  load  for  MIMO  DSP  is  proportional  to  the 
modal  group  delay  [8].  Interestingly,  coupled  multi-core  fibers 
can  be  designed  to  have  reduced  modal  group  delays  and/or 
larger  effective  areas  in  comparison  with  few-mode  fibers  (FMF) 
[9] .  In  a  recent  experiment,  the  modal  delay  was  also  found  to 
depend  sublinearly  on  transmission  distance  in  the  presence 
of  strong  supermode  coupling  [10].  It  is  thus  imperative  to 
study  the  properties  of  optical  supermodes  in  a  systematic  man¬ 
ner.  Such  study  can  help  predict  the  mode  behavior  in  various 
coupled  multi-core  waveguides,  including  single-mode  or  few¬ 
mode  coupled  multi-core  fibers  [10],  [11],  directional  couplers 
or  the  slow-tapering  regime  of  the  multi-core  lanterns  [12],  [13]. 

In  this  paper,  we  analyze  and  characterize  the  modes  of  multi¬ 
core  structures  that  could  be  potentially  useful  in  MDM  systems. 
For  structures  with  single  mode  constituent  elements,  we  in¬ 
vestigate,  in  Section  II,  supermode  arrangements  in  multi-ring 
honeycomb  lattices.  Supermode  structures  with  multimode  con¬ 
stituent  elements  were  first  introduced  in  [14].  In  Section  III, 
we  provide  analytical  description  for  higher-order  supermodes 
in  such  array  configurations,  including  linear,  ring  and  square 
array  lattices.  Finally  in  Section  IV,  we  present  an  experimental 
observation  of  higher-order  supermodes  in  a  3 -core  fiber  array. 

II.  Supermodes  in  Multi-Ring  Honeycomb  Structures 

Honeycomb  waveguide  arrangements  have  been  intensely  in¬ 
vestigated  over  the  years  [15]-[22].  Their  close-packed  geom¬ 
etry  makes  them  ideal  for  multi-core  fiber  application.  In  this 
section,  we  develop  a  general  methodology  capable  of  analyz¬ 
ing  supermodes  in  multi-ring  honeycomb  waveguide  lattices, 
consisting  of  any  number  of  layers.  In  addition  to  hexagonal 
systems,  this  versatile  technique  can  be  used  in  other  regu¬ 
lar  polygonal  configurations.  In  brief,  our  approach  is  based 
on  separately  formulating  an  eigenvalue  problem  for  in-phase 
supermodes  and  another  for  modes  with  orbital  angular  mo¬ 
mentum.  This  can  be  achieved  by  appropriately  grouping  the 
constituent  waveguide  elements  based  on  their  topology.  Here, 
this  approach  is  explicitly  demonstrated  in  honeycomb  arrange¬ 
ments  involving  7,  19,  and  37  sites. 

We  begin  our  analysis  by  assuming  that  in  all  cases,  each 
waveguide  element  is  assumed  to  be  cylindrical  (of  radius  a) 
and  single-moded,  i.e.,  it  only  supports  the  LP0 1  mode,  with 
a  propagation  constant  .  The  coupling  coefficient  k  between 
waveguide  channels  can  then  be  obtained  using  coupled-mode 
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In  Eq.  (2.2)  we  only  account  for  nearest  neighbor  interactions 
in  this  hexagonal  system.  By  introducing  a  gauge  transformation 
Un  =  un  exp  (ifioz)  and  by  adopting  a  dimensionless  coordi¬ 
nate  Z  =  kz  one  arrives  at 
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+  uq  +  {un+ 1  +  'Un-i)  —  0-  (2.3b) 

In  all  cases,  the  modes  can  be  re-expressed  in  terms  of  ac¬ 
tual  quantities  and  coordinates  in  a  straightforward  fashion.  In 
solving  this  problem,  we  look  for  discrete  supermodes,  hence¬ 
forth  denoted  as  DSim ,  where  the  discrete  index  l  represents 
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Fig.  1.  (a)  A  7-core  hexagonal  lattice.  All  single  mode  waveguide  elements 

interact  with  their  nearest  neighbors  with  a  coupling  strengths,  (b)  Real  part  of 
the  held  distribution cos  (Qn)) associated  with  the  DS\\ supermode. 


theory.  If  the  distance  between  core  centers,  i.e.,  the  core  pitch 
is  D ,  this  coupling  strength  is  given  by  [23],  [24]: 


(a) 


(b) 


where  A  =  (nc  —  ns)  /nc  is  the  waveguide  index  difference 
and  nc,  ns  are  the  core  and  cladding  refractive  indices,  re¬ 
spectively.  V  =  koancV2A  is  the  dimensionless  V  number  in¬ 
volved  in  the  eigenvalue  problem  describing  the  fundamental 
mode  LPoi  of  the  waveguide,  e.g.,  V  J\  ( U )  /Jo  (U)  =  WK\ 
(W)  /Kq  (W).  In  the  last  equation,  U  =  a  (A^n2  —  P2)1^2  ,W 

=  a  ( P 2  —  ^ons)^2 ’  where  Ki  ( x )  and  J/  (x)  are  Bessel  func¬ 
tions  of  order  /.  The  propagation  constant  P  can  be  determined 
from  the  eigenvalue  problem  by  keeping  in  mind  that  V2  = 

U 2  +  W2. 

In  this  section,  we  will  exemplify  our  method  in  the  case  of  a 
7-,  19-,  and  37-core  hexagonal  lattice.  Even  though  both  the  7- 
and  19-core  systems  have  been  previously  analyzed  using  either 
direct  schemes  or  linear  algebraic  methods  [18],  [21],  [22],  here 
we  will  still  apply  our  methodology  to  these  same  geometries 
for  demonstration  purposes. 

A  7 -core  arrangement  is  shown  in  Fig.  1(a).  In  this  multi-core 
system  the  modal  fields  evolve  according  to: 


Fig.  2.  (a)  A  19-core  hexagonal  lattice.  The  waveguide  elements  are  grouped 

in  four  different  families  (a,6,c,  d)  indicated  with  different  colors,  (b)  By 
excluding  the  central  core,  the  remaining  elements  of  these  families  are  in 
addition  grouped  in  each  hexagonal  sector  in  order  to  obtain  supermodes  having 
angular  momenta. 


an  orbital  angular  momentum  and  m  is  the  root  number  of 
the  corresponding  eigenvalue  equation.  To  obtain  the  seven  su¬ 
permodes  of  the  structure  depicted  in  Fig.  1(a),  we  first  seek 
in-phase  modes  of  the  type  DSom  that  represent  to  some  ex¬ 
tent  discrete  analogues  to  the  continuous  LPom  modes  in  stan¬ 
dard  cylindrical  waveguides  [24,  25].  Under  this  assumption 
and  by  exploiting  the  symmetry  of  this  problem,  one  can  iden¬ 
tify  two  species,  ua,  u &  where  ua  =  uo  and  =  u\  =  — 

—  uq.  From  Eq.  (2.3)  we  find  that  idua/dZ  +  =  0 

and  idujj/dZ  +  2 +  ua  =  0.  To  obtain  the  dispersion  rela¬ 
tion  of  these  in-phase  modes  we  assume  solutions  of  the  type 
ua  =  vPaelXZ  and  =  u^eaz  that  directly  lead  to  an  algebraic 
eigenvalue  equation  for  the  normalized  propagation  constant 
A, A2  —  2 A  —  6  =  0,  with  solutions  k  =  1  d=  y/7.  As  we  will  see 
these  two  eigenvalues  correspond  to  the  two  in-phase  super¬ 
modes  DSoi  and  DS02  •  From  here,  the  eigenvectors  [u®  ]  can 

be  readily  deduced  from  the  corresponding  eigenvalue  problem. 

We  next  consider  supermodes  involving  angular  momentum. 
In  this  respect,  the  modal  fields  in  the  outer  ring  vary  according 


to  un  =  Beaz  e1®71 .  In  addition,  we  assume  that  the  field  in  the 
central  element  is  zero  (uq  =  0).  In  this  case 


.  dun 

1  1“  ^n  +  1  T  ^n  —  1  =  0- 


(2.4) 


Given  that  any  solution  of  this  latter  equation  must  repeat 
itself  after  six  elements,  (e.g.,  6Q  =  then  upon  substitution 

of  un  in  Eq.  (2.4)  we  deduce  that 


X  =  2  cos  (Q) 


(2.5) 


where  Q  =  In/ 3  and  the  angular  momentum  index  takes  values 
from  the  set  l  =  1,  2, . . . ,  5.  This  eigenvalue  relation  (Eq.  (2.5)) 
can  now  provide  the  remaining  five  supermodes  DSn.  Note 
that  for  the  form  assumed,  Yln= 1  un  =  ®  which  is  consistent 
with  the  assumption  that  the  DSim  fields  are  zero  at  the  center 
(uo  =  0)  when  l  is  finite.  Fig.  1(b)  illustrate  the  field  profile  of 
the  .Das'!  1  supermode. 

This  same  principle  can  now  be  applied  to  other  honeycomb 
lattices  having  more  rings.  As  a  next  example  we  consider  a  19- 
element  honeycomb  double  ring  structure  as  shown  in  Fig.  2(a). 
Given  the  topology  of  this  configuration,  the  elements  involved 
are  now  classified  in  four  groups  denoted  as  a,  b ,  c,  andd — 
illustrated  with  different  colors.  As  in  the  previous  section, 
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Fig.  3.  Real  part  of  the  modal  field  distributions  in  a  19-core  lattice,  associated 

with  (a)  the  DSqi  (b)DSo2  (c)  DS 03  (d)  DSqa  supermodes. 


(a)  (b)  (c) 


Fig.  4.  Real  part  of  the  modal  field  distributions  in  a  19-core  lattice, 
associated  with  supermodes  having  angular  momenta  (a)  DS\\  (b)  DS21  and 
(c)  DS31  mode. 


TABLE  I 

The  list  of  Eigenvalues  Obtained  from  the  Cubic  Equation 
for  the  Modes  with  Angular  Momentum 


Q 

Eigenvalues 

l 

1  =  3.354,  -0.476,  -1.877 

2 

X  =  y/3,  1,  -y/S 

3 

X  sf  -1  +  >/2,0,  —1  —  \/2 

4 

x  =  Vs,  i, -Vs 

5 

X  =  3.354,  -0.476,  -1.877 

because  of  the  hexagonal  geometry,  all  waveguides  interact 
only  with  their  nearest  neighbors  with  a  coupling  strength  k. 
Higher  order  interactions  are  here  ignored  since  the  coupling 
coefficient  tends  to  exponentially  decrease  with  distance.  We 
now  look  for  in-phase  supermodes  with  l  =  0.  In  this  situation, 
the  corresponding  evolution  equations  for  the  modal  fields,  for 
these  four  groups,  read  as  follows: 
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Again,  the  in-phase  eigenvalue  problem  can  be  obtained 
using  the  eigenvectors  uaj,c,d  =  u a  b  c  de%XZ  from  where  we 
find  that  A.4  —  15A.2  —  2A.3  +  24  =  0.  The  four  roots  of  this 
algebraic  equation  are  given  by  [A4  =  4.8715,  A. 2  =  1.2267, 
A, 3  =  —1.6215,  A.4  =  —2.4767]  and  correspond  to  the  DSqi, 
DS02,  DSos  and  DS04  in-phase  supermodes.  Having  found  the 
eigenvalues,  the  corresponding  eigenvectors  u®,  v%,  u®,  u°d  can 
be  evaluated  from  the  matrix  of  Eq.  (2.6).  Fig.  3  depicts  the  real 
part  of  modal  field  profile  of  the  DSom  supermodes.  The  re¬ 
maining  15  modes  are  only  possible  if  the  angular  momentum  l 
is  finite.  In  this  case  again  the  field  in  the  central  element  is  zero 
(ua  =  0).  In  view  of  the  b ,  c,  and  d  groups  assigned,  Fig.  2(b) 
suggests  that 


equation  are  tabulated  in  Table  I  for  each  angular  momentum 
1.  In  this  way  the  remaining  15  supermodes  can  be  generated 
through  the  corresponding  eigenvectors  of  the  matrix  of  Eq. 
(2.7).  As  an  example,  the  real  part  of  field  profiles  of  modes 
DSn,  DS21  and  DSsi  are  shown  in  Fig.  4(a),  (b),  and  (c). 

Other  variations  of  these  hexagonal  lattices  can  also  be  in¬ 
vestigated  using  this  method.  For  example  let  us  consider  again 
the  19-site  configuration  when  this  time  the  central  waveguide 
is  now  completely  removed.  Using  the  same  grouping  outlined 
above,  the  in-phase  modes  satisfy 
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and  hence  the  eigenvalue  equation  is  given  by  A.3  —  2A.2  —  9  A  = 
0,  having  three  roots  A.  =  —2.1623,0,4.1623.  One  can  show 
that  the  remaining  15  modes  ( DSim )  of  this  1 8-site  arrangement 
are  identical  to  those  obtained  for  the  19-site  lattice  when  l  = 
1,2,. ..,5. 

As  a  final  example,  we  investigate  a  37  (three  layer) 
hexagonal  lattice  consisting  of  single-mode  waveguides  as 
shown  in  Fig.  5(a).  Because  of  their  topology,  the  elements 
are  now  grouped  in  seven  families  (a,  b ,  c,  d,  e,  /,  and  g)  as 
illustrated  with  different  colors  in  Fig.  5(a).  Assuming  in-phase 
modes,  we  find 


iiihn  +  ubn  +  1  Jrubn_l  +uCn  +  udn  +udn_l  =0  (2.7a) 

iuCn  +ubn  +  udn  +udn^t  =0  (2.7b) 

iiidn  +  uK  +  ubn  +  1  +uCn  +uCn  +  l  =0  (2.7c) 

where  iibn  =dubn/dZ ,  etc.  By  employing  the  ansatz 
[Ubn ,  uCn  ,ud  J  =  [uln ,  u°Cn ,  u°dn  ]  eiQneixz  in  equation  (2.7), 
we  arrive  at  A.3  —  2A2cos  (Q)  —  X  [4  cos  (Q)  +  5]  —  4  sin2 
(Q)  =  0.  Again,  periodicity  demands  that  Q  =  In/ 3  where  l 
is  1,  2, . . . ,  5.  The  three  eigenvalues  obtained  from  this  cubic 
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Fig.  5.  (a)  37-core  hexagonal  lattice.  The  waveguide  elements  are  now 

grouped  in  seven  different  families  (a,  6,  c,  <i,  e,  /,  g)  indicated  with  different 
colors,  (b)  By  excluding  the  central  core,  the  remaining  elements  of  these  fam¬ 
ilies  are  grouped  in  each  hexagonal  sector  in  order  to  obtain  all  the  supermodes 
with  angular  momenta. 

Seven  supermodes  DSom  can  then  be  obtained  from  the  re¬ 
sulting  algebraic  eigenvalue  equation  A7  —  2A.6  —  23A,5  —  2A,4 
+  124A,3  +  164A,2  +  46 A.  —  12  =  0.  The  seven  roots  of  this  lat¬ 
ter  equation  are  given  by  A.  =  5.3492,2.9752,0.1594,  —0.9524, 
—1.0000,  —1.8567,  —2.6748.  As  before,  the  remaining  30  su¬ 
permodes  exhibit  angular  momentum.  These  can  be  determined 
from  the  field  evolution  equations 


Fig.  6.  Real  part  of  the  modal  field  distributions  in  a  37-core  lattice, 
with  supermodes  (a)  DS03  (b)  DSqa  (c)  D S2i  (d)  D S23  (e)  DS32 
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J_ i L 


(2.11) 


where  Q  =  ln/3  with  l  =  1,  2, . . . ,  5.  The  6x5  eigenvalues 
of  these  DSim  supermodes  (carrying  angular  momentum)  can 
be  directly  obtained  from  Eq.  (2.11)  along  with  their  eigenvec¬ 
tors.  The  real  part  of  modal  field  profiles  of  the  DS03,  Z)So4, 
DS21,  DS23,  DS32  and  DS35  is  depicted  in  Fig.  6.  Finally  the 
eigenvalue  distributions  and  degeneracies  of  a  7-,  19-,  37-core 
lattice  are  provided  for  comparison  in  Fig.  7.  We  would  like  to 
note  that  in  all  cases  our  results  (based  on  coupled-mode  the¬ 
ory)  are  in  excellent  agreement  with  finite  element  simulations 


Fig.  7.  Eigenvalue  distribution  diagrams  for  (a)  a  7-core  (b)  19-core  and  (c) 
37-core  hexagonal  lattice. 

based  on  weakly  guiding  structures.  The  method  described  here 
can  be  used  in  a  similar  fashion  to  analyze  any  other  multi-ring 
hexagonal  lattice.  For  example,  in  the  case  of  a  61 -core  hexago¬ 
nal  lattice  (involving  four  layers),  the  grouping  of  the  elements 
will  lead  to  1 1  in-phase  DS^m  supermodes  (whose  eigenvalues 
are  obtained  from  a  1 1th  order  polynomial)  while  10  x  5  DSim 
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Fig .  8 .  Coupling  between  two  arbitrarily  oriented  LP\  \  degenerate  modes  as  a 
function  of  their  initial  angles  tpi  and  (f2  •  The  arrow  points  to  the  position  of  zero 
coupling  k((£i  =0,(^2  =  7t/2)  =  0  between theL P\ \x  andL P\ \y modes;  the 
inset  shows  the  coordinate  system  of  the  two  cores  including  reference  axes  and 
two  initial  angles. 


LPqi  supermodes.  However,  the  higher-order  LPpq  ( p  >  0) 
“core  modes”  can  also  form  supermodes.  Moreover,  these  super¬ 
modes  are  strongly  affected  by  the  geometrical  distribution  of 
the  cores  within  the  MCF,  because  the  coupling  of  LPpq  (p  >  0) 
modes  between  two  cores  varies  significantly  with  initial  angles 
of  the  modes.  We  name  these  angle-dependent  supermodes  as 
“higher-order  supermodes.”  Unlike  the  higher-order  core  modes 
(p  or  p  <  0),  higher-order  supermodes  specifically  refer  to  the 
supermodes  formed  because  of  coupling  among  “core  modes” 
with  p  >  0.  In  this  section,  we  begin  the  study  of  higher-order 
supermodes  from  the  simplest  two-core  structure  by  demon¬ 
strating  how  angle-dependent  coupling  influences  the  eigen¬ 
mode  formation.  Then,  we  extend  the  theory  to  more  complex 
geometrical  structures,  using  symmetries  to  produce  analyti¬ 
cal  formulas  for  higher-order  supermodes.  More  specifically, 
we  derive  formulas  for  commonly  used  formations,  including 
linear-array,  grid-array  and  ring-array  structures. 


supermodes  with  angular  momentum  can  result  from  a  10th 
order  polynomial  involving  Q  =  In/  3 where  l  =  1,2,. ..,5. 


III.  Theory  of  Higher-Order  Supermodes 


When  the  guiding  power  is  increased  for  an  isolated  core, 
it  can  guide  high-order  modes.  For  optical  fibers  with  low  re¬ 
fractive  index  contrast,  the  vector  modes  are  weakly  guided  and 
therefore  can  be  treated  as  linear  polarization  (LP)  modes.  In 
essence,  the  LP  modes  are  scalar  approximations  of  the  vector 
mode  fields  and  contain  only  one  transverse  field  component 
given  by, 


Epq(r,0 )  =  < 


U 


Jp  {npq  •  r ) 
Jp  (fcpq  ’  cl) 

jj  Ep  (7 pq  ’  r ) 
Kp  {K pq  ’  a) 


cos (p0  +  <£>), for r  <  a 


cos  (pO  +  p),  for r  >  a 


(3.1) 


where  U  is  the  complex  amplitude,  J  and  K  are  Bessel  func¬ 
tions  of  the  first  and  second  kind,  a  is  the  core  radius,  p  is  an 
arbitrary  start  angle,  and  p  is  a  non-negative  integer  referred  to 
as  the  azimuthal  mode  order.  For  the  same  p  value,  npq  and  jpq 
can  take  on  multiple  discrete  values  determined  from  the  dis¬ 
persion  relation.  Therefore  it  is  appropriate  to  label  each  of  the 
aforementioned  terms  with  another  non-negative  integer  corre¬ 
sponding  to  the  number  of  times  the  field  crosses  zero  along 
the  radial  direction.  Thus,  the  LP  modes  can  be  labeled  as 
LPpq.  The  orthogonal  field  component  is  the  same  as  Eq.  (3.1) 
therefore  the  results  presented  here  applies  to  either  polariza¬ 
tion.  According  to  Eq.  (3.1),  when  p  >  0,  the  modal  field  has 
an  additional  degree  of  freedom  in  the  azimuthal  direction  be¬ 
cause  the  start  angle  p  is  arbitrary,  which  generates  spatially 
orthogonal  modes  even  if  both  p  and  q  are  the  same.  For  ex¬ 
ample  the  LPnx  and  LP\\y  modes  correspond  to  the  cases  of 
p  =  0  and  p  =  7r/2,  respectively,  shown  in  Fig.  8.  This  is  very 
different  from  the  case  of  p  =  0  where  the  LPqq  modal  field  is 
azimuthally  uniform  and  has  no  spatial  degeneracy. 

To  the  best  of  our  knowledge,  until  now  the  study  of  super¬ 
modes  has  been  limited  to  the  supermodes  composed  of  these 
azimuthally  uniform  LPpq(p  =  0)  “core  modes,”  particularly 


A.  Higher-Order  Supermodes  in  a  Two-Core  Structure 

As  described  above,  the  coupling  between  two  LPpq  (p  >  0) 
modes  strongly  depends  on  the  initial  angles  of  both  modes.  In 
order  to  focus  on  this  angular  dependence,  the  LPpq  (p  >  0) 
modal  field  of  Eq.  (3.1)  is  written  as  a  function  of  the  initial 
angle  p 

E(p )  =  Fr  •  cos  (pO  —  p) 

=  Ex  •  cos  p  +  Ey  •  sin  p  (3.2) 

where  Fr  represents  the  radial  component  of  the  modal  field, 
Ex  =  E(p  =  0)  and  Ey  =  E(p  =  f)  correspond  to  the  modal 
fields  aligned  with  the  horizontal  axis  x  and  vertical  axis  y. 
A  coordinate  system  is  selected  where  the  azimuthal  reference 
axis  is  parallel  to  the  edge  of  the  graph  formed  by  vertices  at 
the  center  of  the  isolated  cores  as  shown  in  the  inset  of  Fig.  8. 
The  coupling  coefficient  between  two  LPpq  (p  >  0)  modal  fields 
with  initial  angles  of  p\  and  p<i ,  respectively,  in  the  two  identical 
cores  1  and  2  is  given  by 

K'(Ei(p1),E2(p2))  =5®  —  /  /  £0  (nCoreB  -  nClad) 

^  J  J core2 

1)  •  E2(<P2)dxdy 

=  Kx  •  COS  pi  COS  P2  +  Ky 

•  sin  pi  sin  p^  +  •  •  •  nxy 

•  cos  pi  sin  P2  +  Kyx  •  sin  pi  cos  p^ 

(3.3) 

where  nx  —  /^(-Ei,#,  E^yX)^  ny  —  K(^Eiy^E2^y),  cLxy  — 
k(Ei,x,  E2:y),  nyx  =  k{Ei^^  E2,x).  The  modal  fields  are 
normalized,  8$  and  uj  are  vacuum  permittivity  and  angular 
frequency,  ncove2  and  nciad  are  the  refractive  indices  of 
core  2  and  the  cladding,  respectively.  Fig.  8  represents  the 
typical  behavior  of  the  coupling  coefficient  as  the  initial 
angles  of  two  LPU  modes  change,  according  to  Eq.  (3.3). 
nx  >  riy  >  0  is  observed  as  a  result  of  the  specific  geometrical 
distribution  of  modal  fields  with  respect  to  the  cores.  More 
importantly,  coupling  vanished  for  two  specific  initial  angles, 
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i.e.,  Kxy  =  Kyx  =  0,  as  illustrated  in  by  the  arrow  in  Fig.  8.  This 
phenomenon  occurs  because  the  mirror-reversal  of  two  modal 
fields  across  an  axis  parallel  to  the  reference  axis  should  have 
the  same  coupling  coefficient,  which  can  be  mathematically 
elaborated  as 


^xy 


nil  £o(nCor  eB  nClad) 

^  J  J  core2 

F*  cos(#i)  •  F2,r  sm(02)dxdy 


=  -tt 

2  J  J Core2 


5o(^Cor  qB  nClad) 


•F*r  cos( — 0\ )  •  F2,r  sm(—02)dxdy 


—  foyx  —  0* 


(3.4) 


With  this  result,  Eq.  (3.3)  can  be  immediately  simplified 
as 


ft(i?l((^l),i?2((^2))  =  kx  •  COS  (fi  cos 

+  Ky  •  sin  ipi  sin  y?2 .  (3.5) 


Similar  to  the  LPqi  supermodes,  higher-order  supermodes 
can  be  investigated  using  coupled-mode  theory  under  the  as¬ 
sumption  of  weak  coupling.  For  the  simplest  case  of  the  two-core 
structure,  the  basis  set  for  each  core  is  strategically  chosen  as 
the  Ex  and  Ey  pair,  even  though  it  could  have  been  any  orthog¬ 
onal  set  with  an  arbitrary  initial  angle.  Then  the  interactions 
between  the  “core  modes”  can  be  described  by  the  following 
coupled-mode  equation  in  a  matrix  form 


i-^-U  +  MU  =  0, 

dz 


(3.6) 


where  U  =  [u\jX  U2yX  u\iV  U2,y]T \um,x  or  um^y  are  the  com¬ 
plex  amplitude  of  the  horizontally  or  vertically  aligned  modal 
field  of  the  mth  core  after  adopting  the  gauge  transforma¬ 
tion  um^x/y  =  Um,x/y  •  elP°z  in  which  /?0  is  the  propagation 
constant  of  any  isolated  core  mode  LPpq ;  and  hence  M  is  a 
4x4  coupled  matrix 

0  Kx  0  l^xy 


M  = 


K,x  0  tZ xy 

0  l^yx  6 


0 

Ky 


(3.7) 


Kyx  0  Ky  0 


The  standard  procedure  of  solving  the  coupled-mode  equa¬ 
tion  is  to  diagonalize  the  coupled  matrix  Q~lMQ  =  A  where 
the  eigenvalues  A  give  the  normalized  propagation  constants  of 
the  supermodes.  The  corresponding  eigenvectors  U'  =  Q  lU 
describes  the  supermode  field  amplitudes  and  the  row  vectors 
of  A  =  Q1  represents  the  amplitude  coefficients  of  superposi¬ 
tion  of  core  modes  in  forming  the  supermodes.  In  the  supermode 
basis,  the  coupled-mode  equation  reduces  to 

Pu'  +  AIJ'  =  0.  (3.8) 

dz 

For  this  particular  coupled-mode  equation,  one  can  apply 
Kxy  =  Kyx  =0  first  and  then  divide  the  equation  into  two 
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Fig.  9.  Modal  fields  of  LPoi  supermodes  (a  and  b)  and  LP\\  supermodes 
(c)-(f)  of  a  basic  two-core  structure.  In-phase  supermodes  are  shown  in  (a),  (c) 
and  (d)  while  out-of  phase  supermodes  are  shown  in  (b),  (e)  and  (f). 
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Fig.  10.  Modal  fields  ofLPoi  supermodes  (a)-(d)  and  LPn  supermodes  (e)- 
(1)  of  a  4-core  linear-array  structure. 


reduced  ones  as  follows: 


.  d 

ldz 

.  d 
lTz 


^l,x 

+ 

"0 

^ X 

^l,x 

=  0 

(3.9a) 

^2,x 

_K>x 

0 

^2,x 

ul,y 

+ 

"0 

Ky 

ul,y 

=  0. 

(3.9b) 

_U2 ,y  _ 

_Ky 

0 

_u2,y_ 

The  above  equations  mean  that  Ex  and  Ey ,  the  horizontally 
and  vertically  aligned  “core  modes”  do  not  “talk  to”  each  other 
and  their  formation  of  higher-order  supermodes  can  be  solved 
separately.  Eqs.  (3.9a)  and  (3.9b)  turn  out  to  be  very  similar  to 
those  for  the  case  of  LP0i  supermodes.  The  amplitude  coeffi¬ 
cients  of  the  higher-order  supermode  corresponding  to  the  row 
vectors  of  A  =  Q  1  =  ^  [1, 1;  1,  —1]  are  the  same  as  those  for 
the  LPqi  supermodes.  For  higher-order  supermodes  composed 
of  Ex  and  Ey  respectively,  as  shown  in  Fig.  9,  the  propagation 
constants  of  the  four  higher-order  supermodes  f3  =  fio  +  A  are 
calculated  to  be  /3q  ±2  kx  and  /3q  ±2  Ky,  respectively. 


B.  Higher- Order  Supermodes  in  Linear-Array  Structures  with 
Translational  Symmetry 

A  slightly  more  complex  structure  is  the  linear  array  where 
the  cores  are  linearly  aligned  with  each  other  with  equal  core- 
to-core  distances  as  shown  in  Fig.  10.  In  this  case,  couplings 
between  non-adjacent  cores  are  expected  to  be  negligible.  Let 
us  start  with  the  LPqi  supermodes,  which  do  not  have  angular 
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dependence.  One  could  solve  the  eigen-problem  of  the  coupled¬ 
mode  equation  as  described  in  the  last  section  and  obtain  the 
supermodes.  The  alternative  method  is  to  find  the  solutions 
that  satisfy  the  boundary  conditions  uq  =  iz/v+i  =0  [24].  As  a 
result,  the  amplitude  of  the  rath  core  within  any  LPoi  supermode 
can  be  described  as 

Am  =a-eimQ  +  b-e~imQ  ,m  =  1, 2,  -  -  -N  (3.10) 


where  a  and  b  are  coefficients;  Q  is  the  common  phase  acquired 
by  shifting  any  one  lattice  due  to  the  translational  symmetry.  The 
boundary  conditions  applied  to  both  sides  of  the  linear  structure 
(m  m  0  and  m  =  N  +  1) 

An+ 1  =  a  ■  ei(N+l)Q  +  b  •  e~i('N+1)Q 

=  Aq  =  a  ■  ei0  +  b  ■  e~i0  =  0,  (3.11) 

yield  a  +  b  =  0  and  Q  =  l  =  1,2, ...  N.  Therefore  the 
LPqi  supermode  can  be  written  as 


(  Trim  \ 

u'lm  =  Ai,m  ■  Uitm  =  sin  I  -  J  ■  U{,m  ,  TO,  l  =  1,  2,  •  •  •  N 

(3.12) 

where  l  corresponds  to  the  order  of  different  supermodes  and  urn 
represents  the  complex  amplitude  of  the  rath  “core  mode.”  For 
the  Zth  LPqi  supermode,  comparing  the  coupled-mode  equation 


.  d 

i  um 
dz 


+  K(Um- 1  +  Um  +  l)  —  0? 


with  the  supermode  equation 


idzU'm  +XlU'm  =  0 


(3.13) 


(3.14) 


where  k  is  the  coupling  coefficient  between  two  adjacent 
LPqi  “core  modes,”  the  propagation  constant  of  the  Zth  LPqi 
supermode  A  =  A)  +  A  can  be  obtained  as 


A  =  A)  +  2/t  •  cos 


7 rZ 


N  +  1 


,Z  =  1,2- • -TV. 


(3.15) 


For  higher-order  supermodes,  because  only  coupling  between 
adjacent  cores  is  considered,  Ex  modes  would  only  couple  to 
themselves  as  would  Ey  modes,  according  to  nxy  =  nyx  —  0. 
Therefore  the  higher-order  supermodes  of  the  linear-array  struc¬ 
ture  can  be  divided  into  the  Ex-  and  Ey -families;  each  of  them 
can  be  solved  independently  using  the  same  relations  as  used 
for  the  LPqi  supermode.  The  propagation  constants  are  attained 
as 


A,ar 


A)  +  2nx  *  cos 


7 rZ 


N  +  l 


,1  =  1,2- • -TV  (3.16a) 


0i, y  =0o+  2kv  ■  cos  ^  i  J  4  =  1, 2  •  •  •  IV  (3.16b) 

respectively.  An  example  of  the  LP\  \  supermodes  as  well  as  the 
LPqi  supermodes  of  a  4-core  linear- array  structure  is  shown  in 
Fig.  10. 


C.  Higher-Order  Supermodes  in  2D  Rectangular- Array 
Structures  with  Direct  and  Diagonal  Coupling  Interactions 

We  then  analyze  the  supermodes  in  two  dimensional  (2D) 
discrete  waveguide  arrays.  These  are  composed  of  M  x  N 


Fig.  11.  (a)  A  rectangular  array  of  waveguides.  The  LPoi  mode  in  each 

waveguide  cross-talks  with  nearest  neighbors  along  the  horizontal  and  vertical 
directions  as  well  as  with  diagonal  elements,  (b)  Coupling  interactions  in  this 
same  array  when  each  element  involves  instead  the  LP\\  mode. 


identical  optical  waveguides  arranged  in  a  rectangular  geometry. 
The  distance  between  cores  is  taken  here  to  be  D.  For  gener¬ 
ality,  we  also  include  in  our  analysis  not  only  nearest  neighbor 
interactions  k  (occurring  along  the  horizontal  and  vertical  di¬ 
rections)  but  also  higher-order  couplings  nc  taking  place  along 
the  two  diagonals  as  shown  in  Fig.  11(a).  For  weakly  guid¬ 
ing  structures,  the  coupling  strengths  can  be  obtained  from  the 
Eq.  (2.1).  We  first  consider  supermodes  derived  from  the  LPqi 
mode  supported  by  each  waveguide  channel  in  a  square  lattice 
(see  Fig.  1 1(a)).  In  such  an  arrangement,  the  modal  fields  evolve 
according  to  [26] : 


i  dz  H-  A)  Um ,  n  H-  tv  {Um  +  l,n  “f  l,n 
P  Um,n+ 1  P  Lfm,n  —  i)  -f  nc  {Um  +  l,n  +  l  H-  ^4n  +  l,n— 1 

EUm— 1,71+1  P  Lfm— i,n— 1)  =  0  (3.17) 


where  Z7n?m  represents  the  modal  field  amplitude  at  site  n,  m 
in  this  rectangular  array  and  the  discrete  site  indices  take  values 
from  the  sets  n  m  1,  2, . . . ,  N  and  m  =  1,2,...,  M.  To  identify 
the  eigenmodes  of  this  system,  we  look  for  solutions  of  the  type 
Um,n  —  dLd°zurrirn  that  also  satisfy  the  boundary  conditions 
Um, o  =  um,N+ 1  =  UQ:n  =  uM+i,n  =  0  [27].  Based  on  these 
requirements,  one  can  directly  show  that  the  Uki1  supermodes  of 
this  structure  are  U^ln  =  A^i1  nel^QZ elXk^lZ  with  the  amplitudes 
given  by 

4A  =  sin  sin  (3.18) 


where  k  =  1,  2, . . . ,  M  and  l  =  1, . . . ,  N.  In  addition,  the  eigen¬ 
value  associated  with  Uk’1  supermode  is  given  by 


Xk,i  =  2n 


cos 


kir 


cos 


M  +  l 

kir 

M  +  l 


+  cos 


hr 


N  +  l 
hr 


N  +  l 


(3.19) 


Evidently,  altogether  this  array  supports  M  x  N  supermodes. 
Similarly,  one  can  investigate  the  eigenmodes  of  this  same  array 
arising  from  the  LP\\  mode  of  each  waveguide  (see  Fig.  1 1(b)). 
In  a  weakly  coupled  array,  the  LPu  mode  tends  to  orient  it¬ 
self  either  along  the  x  or  y  direction  ( LPllx ,  LPlly).  As  a 
result  the  coupling  strengths  nx  and  ny  are  different  because  of 
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a  i  kir  \  f  In 

+4kc  cos  (  —  „  )  cos 


M  +  l 


N+  1 


where  Am,  A+ i,m  and  A,™+i  are  Kronecker  deltas;  l  and  m 
are  integer  numbers  mod  N.  This  matrix  can  be  diagonalized 
by  a  lattice  Fourier  transform  as  A  =  Q~lMQ  [5].  Q  1  is  a 
unitary  matrix  with  elements  given  as 


(<5 1)im  _  we* 


(3.23) 


which  represents  the  amplitude  of  mth  core  for  the  Zth  LP0i 
supermode.  These  discrete  helical  phases  give  supermodes  an 
appearance  similar  to  optical  vortices  but  in  a  discrete  form 
and  thus  carry  orbital  angular  momentum  [4].  The  propagation 
constant  of  the  Ith  LPqi  supermode  equals  the  corresponding 
eigenvalue  and  can  be  calculated  as 


A  s*  A)  +  2k  •  cos 


■N. 


(3.24) 


Fig.  12.  (a)  Scheme  of  mode  basis  selection  for  LP\  i  supermode  analysis  of  a 

4-core  grid-array  structure;  modal  fields  of  LPq i  supermodes  (a)-(d)  and  LP\\ 
supermodes  (e)-(l)  of  a  4-core  grid-array  structure  computed  using  COMSOL. 


their  respective  overlap  integrals.  Hence,  the  field  evolution  is 
described  by 

i  ~dz  ^  A)^m,n  T-  Kx  +  H-  Urri—l.n) 

~\~Ky  {Um,n-\- 1  H-  — l)  H- 

«c  (^4n  +  l,n-|-l  H-  ^m  +  l,n  —  1  H-  Um  —  l,n  +  l  H-  Um—\  ,n—l )  =  0. 

(3.20) 

The  supermodes  Uk’1  of  this  latter  equation  are  exactly  iden¬ 
tical  in  form  with  those  provided  by  Eq.  (3.18).  In  this  case 
however,  the  corresponding  eigenvalues  are  given  by 

f  kn  \  (  hr 

kb  i  =  2/w  cos  — -  +  2 kv  cos  — - 

\M  +  lJ  y  \N  + 1 


According  to  Eq.  (3.24),  it  is  obvious  that  the  Zth  and 
(N  —  Z)th  supermodes  are  degenerate  modes.  This  has  a  phys¬ 
ical  explanation  from  the  equivalence  between  clockwise  and 
counter-clockwise  mode  orders.  A  unitary  transformation  can 
project  the  vortex-like  basis  of  the  Zth  and  ( N  —  Z)th  super¬ 
modes  into  another  orthogonal  basis.  In  particular,  the  degen¬ 
erate  supermodes  can  be  transformed  into  a  basis  in  real  fields 
with  amplitude  transformation  as  below 


(3.25) 


Bl,m 

1 

"1 

1  " 

_  Bn  —i,m  _ 

=  7i 

1 

-1 

_AN-l,m  _ 

where  B,m  =  A/|r  cos  (^r  ■  m)  and  BN_ 


(3.21) 


In  all  cases  our  results  (based  on  coupled-mode  theory)  are 
in  excellent  agreement  with  finite  element  simulations  using 
COMSOL  as  shown  in  Fig.  12. 

D.  Higher-Order  Supermodes  in  Ring- Array  Structures  with 
Rotational  Symmetry 

Ring-array  structure  is  another  interesting  geometry  for 
higher-order  supermodes  as  it  possesses  rotational  symmetry. 
Here  also,  we  only  consider  coupling  between  adjacent  cores. 
We  start  from  PPoi  supermodes  of  N-core  ring-array  structure. 
In  this  case,  the  coupled  matrix  M  is  a  symmetric  circulant  ma¬ 
trix.  If  A)  and  n  are  defined  as  the  propagation  constant  of  the 
“core  mode”  and  the  coupling  coefficient  between  two  adjacent 
PPoi  “core  modes”  respectively,  the  elements  of  the  N  x  N 
coupled  matrix  are  given  as 

(AT)/m  =  A) Am  +  P*  (A+l,m  H"  A,m  + 1)  (3.22) 


=  y  ]v  cos  nr  •  m)  anu  nN-\,m  =  y  sin 

(¥•-)• 

The  higher-order  supermodes  of  ring-array  structure  are  par¬ 
ticularly  interesting  because  they  can  no  longer  be  divided  into 
Ex-  and  Ey -basis  for  independent  analysis.  Instead,  the  tan¬ 
gential  and  normal  mode  basis  Etm  and  EUm  (m  refers  to  the 
core  number)  are  selected,  which  correspond  to  “core  modes” 
aligned  horizontally  or  vertically  with  respect  to  the  reference 
axis  pointing  from  the  ring  center  to  the  mth  core.  Fig.  13(a) 
shows  such  a  mode  basis  for  a  3 -core  ring-array  structure.  In  that 
case,  2  =  k 2,3  =  ^3,1  is  obtained  because  of  the  rotational 
symmetry.  In  addition,  the  coupling  coefficients  between  the 


LPpq  modes  are  real  so  that 


Kjp.  Therefore  the  number 


of  coupling  coefficients  for  the  N-core  ring-array  structure  is 
reduced  from  2N  to  4  given  by 


-K*sin2(v) 

+  ny  coe 

'2(v)- 

(3.26a) 

,n2  = 

sin  (A)  cos  (A) 

—  Ky  sin  ( -^  ) 

™  (v)  • 

(3.26b) 

= 

kx  sin  (-^)  cos 

-*,sm(A)co8(»), 

(3.26c) 
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ring-array  structure 


Fig.  13.  (a)  Scheme  of  mode  basis  selection  for  LP\  \  supermode  analysis  of  a 
3-core  ring-array  structure;  modal  fields  of  LPq  i  supermodes  (b)-(d)  and  LP\  i 
supermodes  (e)-(j)  of  a  3-core  ring-array  structure  computed  using  COMSOL. 


Then  the  coupled-mode  equation  for  the  N-core  ring-array 
can  be  written  as 


dz 


'Ut  ' 

'Ut  ■ 

Un_ 

+  M 

un_ 

=  0, 


(3.27) 


where  Ut  and  Un  are  column  vectors  of  dimension  N.  the 
coupled  matrix 


M  = 


Mtt  Mtn 

Mnt  Mnn 


Tunable  Laser 
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f  Computer  1 
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Control  J 

Fig.  14.  (a)  Schematic  of  spatially-and-spectrally-resolved  imaging  (S2  imag¬ 

ing)  setup  using  a  tunable  laser  and  a  CCD  camera;  (b)  wavelength-scanning 
results  of  the  multi-path  interference  pattern  for  one  pixel  of  the  CCD  camera 
and  (c)  the  result  of  intensity  vs.  differential  group  delay  (DGD)  after  taking 
Fourier  transform  of  (b)  for  the  same  pixel.  Inset  of  (c)  shows  the  resolved  LP 
modes  after  picking  the  information  for  every  pixel  at  the  corresponding  DGDs 
of  different  modes  and  mapping  them  together. 


where  Sim,  ^+i,m  and  are  the  Kronecker  deltas;  l 

and  m  are  integer  numbers  mod  N.  Diagonalization  of  the 
coupled  matrix  M  is  the  key  to  solving  the  higher-order 
supermodes.  It  takes  two  steps  as  follows  as  described  in 
appendix  A. 

Both  the  eigenvalues  and  eigenvectors  can  be  obtained  after 
successful  diagonalization  of  the  coupled  matrix.  After  plugging 
the  coupling  coefficients  into  the  eigenvalues  and  eigenvectors 
using  Eq.  (3.22),  the  higher-order  supermodes  can  be  divided 
into  (±)  groups  and  solved  as  Eq.  (3.29),  shown  at  the  bottom 
of  the  page,  where 

7^  =  —i  (rj  ±  y/l  +  rf^j 

and 

_  ( Kx  -  Ky)  cot  (ffi) 

V  («x  +  Ky)  sin  (^) 

The  propagation  constants  of  the  higher-order  supermodes  are 
the  eigenvalues 


is  a  rank-2N  Hermitian  matrix  and  each  sub-matrix  is  a  circulant 
matrix  with  elements  described  as 

=  Ktlft2  (^+i,m  +  8i,m+i)  5  (3.28a) 
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{Mnn)lm  =  Kni,ri2  '  + 1  “b  ^Z+l,m)  •  (3.28d) 
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where  l  =  1,  2  •  •  •  N.  Same  as  the  LP0i  supermodes,  the  Zth 
and  (N  —  Z)th  higher-order  supermodes  for  both  (d=)  groups 
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Fig.  15.  (a)  Cross-sectional  view  of  the  fabricated  coupled-3-core  fiber  and 

(b)-(j)  resolved  supermodes  of  the  fabricated  coupled-3-core  fiber  shown  in 
amplitude  (xi)  and  phase  (X2):  (b)  the  fundamental  LPoi  supermode;  (c),  (d) 
the  degenerate  pair  of  LPq  i  supermodes;  (d)  the  fundamental  LP\  \  supermode; 
(e),  (f)  the  first  degenerate  pair  of  LP\\  supermodes;  (g)  the  fourth  LP\\ 
supermode;  (h),  (i)  the  second  degenerate  pair  of  LP\\  supermodes. 

are  degenerate.  Real-field  basis  can  be  obtained  by  applying  a 
unitary  transformation  on  the  current  basis.  An  example  of  the 
LP\  i  supermodes  as  well  as  the  LP0i  -supermodes  of  a  3-core 
ring-array  structure  is  shown  in  Fig.  13(b)— (j). 

IV.  Observation  of  Higher-Order  Supermodes 

A  coupled  few-mode  3 -core  fiber  was  fabricated  in  order 
to  observe  higher-order  supermodes.  The  cross-section  of  the 
fiber  is  shown  in  Fig.  15(a),  where  the  Ge-doped  cores  are 
hexagonally-shaped  as  a  result  of  the  stack-and-draw  fabrication 
process.  The  cladding  diameter  is  ^120  fi m,  and  the  index 
difference  between  the  core  and  cladding  is  ^0.6%.  The  core 
diameter  was  selected  to  be  ^9  fim  to  ensure  that  the  cut-off 
wavelength  of  the  fundamental  modes  is  well  above  1.5  fi m. 
In  addition,  the  core  pitch  was  chosen  to  be  ^1 1.5  /jm,  a  small 
enough  distance  to  allow  for  strong  coupling.  Therefore  the  fiber 
supports  both  LPoi  and  PPn  supermodes  around  1.5  fim. 

A  spectrally-and- spatially  resolved  imaging  (S'2  imaging) 
setup  was  built  to  acquire  the  modal  fields  from  the  fiber.  The 
S2  imaging  procedure,  first  developed  by  Nicholson  et  al.  [28], 
is  a  technique  specially  designed  for  quantifying  mode  content 
in  fibers.  The  principle  is  to  spatially  resolve,  in  a  point-by- 
point  fashion,  the  spectral  multi-path  interference  patterns  pro¬ 
duced  by  mode  beatings.  The  imaging  setup  can  either  involve 
a  broadband  source  accompanied  by  a  spatial- scanning  system 
or  a  tunable  laser  source  accompanied  by  a  CCD  camera.  For 
the  purposes  of  this  experiment,  we  chose  the  latter  setup  as 
shown  in  Fig.  14(a)  due  to  its  fast  scanning  speed  and  high  fre¬ 
quency  resolution.  Both  the  tunable  laser  and  the  camera  scan  in 
sub-seconds  and  can  be  easily  synchronized  by  computer  con¬ 
trol.  In  addition,  the  tunable  laser  supports  a  fine  resolution  of 
0.0001  nm  with  a  tuning  range  of  1510-1640  nm.  Light  coupled 
from  the  tunable  laser  source  was  launched  into  approximately 
7  m  of  the  fabricated  3 -core  fiber  through  free  space.  The  out¬ 
put  of  the  3 -core  fiber  is  imaged  onto  the  camera  through  a  4-f 


optical  system.  A  polarizer  was  added  in  front  of  the  camera  to 
guarantee  a  single  polarization.  During  the  scanning  operation, 
the  laser  wavelength  was  incremented  in  discrete  steps  and  a 
mode-interference  image  was  captured  by  the  camera  at  each 
step.  The  measured  S 2  results  can  be  expressed  as  wavelength- 
interference  patterns  for  every  image  pixel.  An  arbitrary  exci¬ 
tation  can  be  decomposed  into  different  fiber  modes,  including 
the  lossy  modes  which  attenuate  quickly  after  a  few  meters  or 
gentle  bends.  Because  modes  travel  at  different  group  velocities, 
if  three  guided  modes  are  assumed  to  be  excited  with  amplitudes 
A,  B,  C  and  phase  <pa,  fo,  <j)c,  the  measured  intensity  at  each 
pixel  location  (x,  y )  would  be 

I{w)  \(x,y)  =  I 

+  B(x,y)  ■  eiMx'y) 

+  C(x,y)  ■  eiMx’v)  .e-iw-T“\ 

=  (A2  +  B2  +  C2)  +  2AB 

•COs(A0oi)  -LO  -Tab) 

+  2 AC  ■  COs(A 4>ac  -  u  ■  Tac) 

+  2BC  •  cos(A(/>6c  -  u>  ■  The),  (4-1) 

where  c o  is  the  angular  frequency  and  represents  the  differ¬ 
ential  group  delay  between  mode  i  and  mode  j.  An  example  of 
a  measurement  at  a  particular  pixel  location  is  plotted  as  shown 
in  Fig.  14(b).  Using  that  data,  one  can  extract  all  the  mode 
information  by  simply  applying  a  Fourier  transform  as  follows: 

i(j)  \(x,y)  =  ( A 2  +  B2  +  C 2)  •  5(t)  +  ... 

AB  •  [elA(/)ab  •  8(t  +  rab )  +  e~lA^ab  •  8(r  -  rah)\ 

+  AC  •  [elA(t>ac  •  8(t  +  Tac)  +  e~lA(t)ac  •  S(r  -  rac)] 
+  BC  •  [eiA+>°  ^(T  +  r6c)+e-lA^  A(r-rbc)\  . 

(4.2) 

Fig.  14(c)  represents  the  same  pixel  location  after  Fourier 
transform,  where  every  peak  corresponds  to  each  mode  beating. 
Therefore,  the  amplitude,  phase,  and  group  delay  of  the  modes 
can  be  resolved  from  the  peak  strength,  phase  and  locations. 
Using  this  method,  each  mode  profile  can  be  extracted  from  the 
series  of  images  by  performing  the  Fourier  transform  on  each 
pixel  as  shown  in  the  inset  of  Fig.  14(c). 

All  the  supermodes  of  the  3 -core  fiber  were  obtained  using 
the  S 2  imaging  method.  One  difficulty  of  the  experiment  is  that 
the  mode  identification  became  very  complicated  due  to  the 
large  number  of  fiber  modes.  Since  every  peak  is  the  result  of 
the  beating  of  any  two  modes,  the  number  of  the  mode-beating 
peaks  scales  quadratically  with  the  number  of  modes.  In  this 
experiment,  the  problem  was  solved  by  using  intentional  off¬ 
set  launching  to  excite  only  one  dominant  mode  with  a  few 
other  modes  for  each  wavelength  sweep.  Then,  multiple  excita¬ 
tions  were  required  in  order  to  resolve  different  modes.  Notice 
that  the  offset  excitation  method  is  selected  here  for  simplicity, 
while  lower  loss  or  more  accurate  excitation  can  be  achieved  by 
using  photonic  lanterns  [13],  [29]  or  spatial  light  modulators. 
S 2  imaging  also  suffers  from  the  inherent  inability  to  separate 
degenerate  modes,  which  share  same  group  delays.  However, 
due  to  the  hexagonal  shape  of  the  cores  in  this  3 -core  fiber, 
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the  degenerate  modes  actually  have  a  subtle  difference  in  group 
delay,  making  it  possible  to  reconstruct  each  of  them  separately. 
Therefore,  a  total  of  nine  modes  supported  by  the  3-core  fiber 
were  successfully  reconstructed  in  both  amplitude  and  phase, 
as  shown  in  Fig.  15.  The  first  three  images  are  of  LPoi  su¬ 
permodes,  in  which  the  field  is  slightly  better  resolved  than  in 
those  of  the  LPn  supermodes.  This  is  because  significantly  less 
power  is  coupled  into  the  higher-order  modes  during  center  or 
offset  excitation.  The  degenerate  higher-order  mode  images  are 
slightly  blurry  because  their  group  delay  difference  is  too  small 
to  allow  them  to  be  clearly  differentiated.  Nevertheless,  the  re¬ 
solved  fiber  modes  are  in  good  agreement  with  the  simulations 
shown  in  Fig.  15. 


V.  Conclusion 

We  have  investigated  supermodes  in  different  waveguide  ar¬ 
ray  structures  based  on  the  weakly-coupled  assumption.  For 
coupled  single-mode  multi-core  waveguides,  supermodes  in 
multi-ring  honeycomb  arrangements,  which  provide  the  highest 
mode  density  per  unit  cross-sectional  area,  have  been  obtained. 
For  coupled  few-mode  multi-core  waveguides,  higher-order  su¬ 
permodes  have  been  analyzed  for  the  linear  arrays,  square  lat¬ 
tices  and  ring  arrays.  General  solutions  of  LPpq  supermodes 
in  those  structures  of  arbitrary  sizes  are  provided.  In  addition, 
an  experimental  observation  of  higher-order  supermodes  has 
been  achieved  for  the  first  time  in  a  coupled  few-mode  3- 
core  fiber.  This  study  enriches  the  concept  of  supermodes  in 
coupled  multi-core  waveguides,  which  may  have  potential  ap¬ 
plications  not  only  in  MDM  systems  [9,  11,  30],  but  also  in 
other  areas  related  to  waveguide  optics,  such  as  optical  phased 
arrays  [31,  32],  beam  combining  [33,  34]  and  fiber  imaging 
systems  [35,  36]. 


Appendix  A 

In  order  to  solve  the  higher-order  supermodes  for  ring-array 
structure,  diagonalization  of  the  coupled  matrix  M  is  the  key.  It 
takes  two  steps  as  follows.  Because  each  sub-matrix  of  M  has 
the  similar  form  as  the  coupled  matrix  of  LP0 1  supermodes,  the 
first  step  is  to  utilize  the  previous  solution  and  diagonalize  the 
sub-matrixes  all  together.  In  order  to  do  so,  a  rank-2N  modal 
matrix  is  defined  as 


=  2  a 


'tl  ,t 2 


(At)M 

(An)i,;  =  2 intun2  ‘  sin 


where  all  the  diagonal  elements  can  be  calculated 
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D  is  a  block  matrix  composed  of  four  diagonal  N  x  N 
sub-matrixes,  which  are  commutable  with  each  other.  Notice 
that  if  (Dtn)i=iQ  =  (Dnt)i=io  =  0,  then  (Dtt)i=i0  =  2a^^2 
cos  (^)  and  ( Dnn)l=lo  =  2 Knii„2  •  cos  (^)are  eigenval- 
ues  already.  In  addition,  a  function  of  a  diagonal  matrix  is  still 
a  diagonal  matrix,  whose  elements  are  the  same  function  of  the 
original  matrix  elements,  i.e.,  f  (D)lm  =  f  ((D)ml)[ 5].  With 
these  properties,  the  second  step  of  diagonalization  can  be  op¬ 
erated  on  D  similar  to  a  simple  2x2  matrix  of  a  ^ 
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The  eigenvalue  matrix  is  solved  as 
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and  the  modal  matrix  is 
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in  which  Dr  =  \f(Du  -  Dnn)2  +  4Di„D„i  when  (. Dtn)l  l 
7^  0.  Therefore  the  higher-order  supermodes  are  the  row  vec¬ 


tors  of  the  matrix  R  1 P  1  where  P  1  =  P*  = 
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